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Abstract 

We report the observation of 823±57 B° and 970±65 B + decays to doubly charmed final states 
D^'D^K, where D [ *> and £> w are fully reconstructed and K is either a K or a K°. We use a 
sample of 82.3 ± 0.9 million BB events collected between 1999 and 2002 with the BABAR detector 
at the PEP-II storage ring at the Stanford Linear Accelerator Center. The 22 possible B decays 
to D { * ] D^K are reconstructed exclusively and the corresponding branching fractions or limits are 
determined. The branching fractions of the B° and of the B + to D^K are found to be 

B(B° -> D^D^K) = (4.3 ± 0.3(stat) ± 0.6(syst)) %, 
B(B+ -> D ( * } D^K) = (3.5 ± 0.3(stat) ± 0.5(syst)) %. 

A search for decays to orbitally excited D s states, B — > Dfj (Dfj — > D^°K + ) is also per- 
formed. 



PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 
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I. INTRODUCTION 



The inconsistency between the measured b — > ccs rate 
and the rate of semileptonic B decays has been a long- 
standing problem in B physics. Until 1994, it was be- 
lieved that the b — > ccs transition was dominated by de- 
cays B — ► D S X, with some smaller contributions from 
decays to charmonium states and to charmed strange 
baryons. Therefore, the b — ► ccs branching fraction was 
computed from the inclusive B — > D s X , B — > (cc)X, 
and B S C X branching fractions, leading to B(b — > 
ccs) = (15.8 ±2.8)% [1]. Theoretical calculations are un- 
able to simultaneously describe this low branching frac- 
tion and the semileptonic branching fraction of the B 
meson [2]. 

As a possible explanation of this problem, it has been 
conjectured [3] that B(b — > ccs) is larger and that de- 
cays of the type B -> D^D^K (X) (where £>W can be 
either a D°, £>*°, D + , or could contribute signifi- 

cantly to the decay rate. This might also include possible 
decays to orbitally-excited D s mesons, B — > D^D s j, 
followed by D s j — > D^K. Experimental evidence in 
support of this picture has been published in the past 
few years. This evidence includes the measured branch- 
ing fraction for wrong-sign D production, averaged over 
charged and neutral B mesons, by CLEO [4] [B(B — > 
DX) = (7.9 ± 2.2)%], and the observation_of a small 
number of fully reconstructed decays B — > D^D^K, 
both by CLEO [5] and ALEPH [6]. More recently, BABAR 
[7] and Belle [8] have reported some preliminary results 
on the evidence for transitions B° — ► D*^ K + with 
much larger data sets. 

B -> decays can proceed through two 

different amplitudes: external IF-emission amplitudes 
and internal W-emission amplitudes (also called color- 
suppressed amplitudes). Some decay modes proceed 
purely through one of these amplitudes while others can 
proceed through both. Figure 1 shows the possible types 
for B -> D^D^K decays. In BABAR, the large data 
sets now available allow comprehensive investigations of 
these transitions. In this paper, we present measure- 
ments of or limits on the branching fractions for all the 
possible B -> D^D^K and B -> D^D^K+ decay 
modes, using events in which both D mesons are fully 
reconstructed. Charge conjugate reactions are assumed 
throughout this paper and branching fractions are aver- 
aged accordingly. 
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FIG. 1: Top: internal IF-emission diagram for the decays 
B — > D^D^K. Bottom: external VU-emission diagram for 
the decays B -> D^D^K. 



II. THE BABAR DETECTOR AND DATASET 

The study reported here uses 75.9 fb _1 of data col- 
lected at the T(4S) resonance with the BABAR detector 
at the PEP-II asymmetric-energy B factory, correspond- 
ing to (82.3 ± 0.9) x 10 6 BB pairs. 

The BABAR detector is a large-acceptance solenoidal 
spectrometer (1.5 T) described in detail elsewhere [9]. 
The analysis described below makes use of charged track 
and 7T° reconstruction and charged particle identification. 
Charged particle trajectories are measured by a 5-layer 
double-sided silicon vertex tracker (SVT) and a 40-layer 
drift chamber (DCH) , which also provide ionization mea- 
surements (dE/dx) used for particle identification. For 
charged tracks with p > 1 GcV/ c, the measured trans- 
verse momentum with respect to the beam axis (px) has 
a resolution a VT such that 

^ = 0.13%p T + 0.45%, (1) 

PT 

where pr is measured in GeV/c. 

Photons and electrons are measured in the electro- 
magnetic calorimeter (EMC). The EMC consists of 
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6580 Thallium-doped Csl crystals, which instrument the 
barrel and forward endcap; the crystals are arranged 
in a quasi-projective geometry. The electromagnetic 
calorimeter resolution ue can be expressed as 

f = *™ © 1.9%, (2) 

where the energy E is measured in GeV. 

Charged particle identification is provided by the av- 
erage energy loss (dE/dx) in the tracking devices and 
by an internally reflecting ring-imaging Chcrcnkov de- 
tector (DIRC). The DIRC comprises 144 quartz bars, 
divided into 12 sectors, which transport the Cherenkov 
light to a water-filled expansion volume equipped with 
10751 photomultiplier tubes. A K/ir separation better 
than four standard deviations is achieved for momenta 
below 3 GeV/c. 

III. B CANDIDATE SELECTION 

The B° and B + mesons are reconstructed in a 

sample of hadronic events for all the possible DDK 

modes, namely B° -> D^-D ( *^K+, D^~D^+K°, 
D (*)o D (*)v K o and B + ^ £(*)0£(* )+ir ; D -(*)o D (*)o K +^ 

D^-D^ + K+. K° mesons are reconstructed only from 
the decays K® — > n + n~. To eliminate the background 
from continuum e+e~ — » qq events, we require that the 
ratio of the second to zeroth Fox- Wolfram moments of 
the event [10] be less than 0.45. 

The K® candidates are reconstructed from two op- 
positely charged tracks consistent with coming from a 
common vertex and having an invariant mass within 
±9MeV/c 2 of the nominal K® mass. For most of the 
channels involving a K®, we require that the K® ver- 
tex be displaced from the interaction point for the event 
by at least 0.2 cm in the plane transverse to the beam 
axis direction. The n° candidates are reconstructed 
from pairs of photons, each with energy greater than 
30 MeV, which are required to have an invariant mass 
115 < m 77 < 150MeV/c 2 . The tt° from a D*° decay 
must have momentum between 70 MeV/c and 450 MeV/c 
in the T(4S) frame, while the 7r° from D° — > K~ir + ir° 
must have energy greater than 200 MeV in the laboratory 
frame. 

The D* candidates are reconstructed in the decay 
modes D*+ -> D°tt+, D*+ -» D+tt°, L>*° -» L>°7r°, and 
D*° — ► D°j. A ±3cr interval around the nominal mass 
difference Am = m(D*) - m(D°) is used to select D* 
mesons, where a is the measured mass difference reso- 
lution and is equal to lMeV/c 2 for D* + — > D°ir + and 
D*+ -» D+tt decays, 1.3MeV/c 2 for D*° -» L>°7r° de- 
cays, and 3.3 MeV/c 2 for D*° -> D° 7 decays. The mode 
D* + — > D + ir° is used only in the reconstruction of decays 
B° -> D*-D*+K° s and B+ -> D*~D*+K+. 

The D° and D + mesons are reconstructed in the de- 
cay modes D° — ► K~n + , K~n + ir°, K~Tr + n~ir + , and 
D + — > i^ _ 7r + 7r + , by selecting track combinations with 



invariant mass within ±2<r of the average measured D 
mass. The average D mass and the D mass resolution a 
used in this selection are fitted from the data itself, using 
a large inclusive sample of D decays. The resolution is 
equal to 7MeV/c 2 for D° -» K~ tt+ decays, 13MeV/c 2 
for D° -> if" 7T+ 7T° decays, 5.7MeV/c 2 for D° -» 
7T+ 7r~ 7T+ decays, and 5.5McV/c 2 for D+ — > if~7r + 7r + 
decays. For modes involving two D° mesons, at least one 
of them is required to decay to K~ir + , except for the 
decay modes D*-D*+K°, D*- D*+ K+ , and D*-D°K+, 
which have lower background and for which all combi- 
nations are accepted. All K and ir tracks are required 
to be well reconstructed in the tracking detectors and to 
originate from a common vertex. Charged kaon identi- 
fication, based on the measured Cherenkov angle in the 
DIRC and the dE / dx measurements in the drift chamber 
and the vertex tracker, is used for most D decay modes, 
as well as for the K + from the B meson decay. 

B candidates are reconstructed by combining one 
one £>(*) and one K candidate. A mass-constrained kine- 
matic fit is applied to all intermediate particles (D*°, 
D* + , D°, D + , K°, 7T°). Since the B mesons are produced 
via e + e~ — > T(4S) — » BB, the energy of the B meson in 
the T(AS) rest frame is given by the beam energy in the 
center-of-mass frame, which is known much more 

precisely than the energy of the B candidate. Therefore, 
to isolate the B meson signal, we use two kinematic vari- 
ables: the difference between the reconstructed energy 
of the B candidate and the beam energy in the center 
of mass frame (AE), and the beam energy substituted 
mass (toes), defined as 

mES = \/(^r) 2 -^ (3) 

where p* B is the momentum of the reconstructed B in 
the T(45) frame. Signal events have toes close to the 
nominal B meson mass, 5.279 GeV/c 2 , and AE close to 
MeV. Due to imperfect modeling of the charged K en- 
ergy loss in the detector material, the central value of 
AE is slightly shifted away from OMeV by an amount 
A-E s hift = (— 5 ± 1) MeV, which is fitted from the data 
themselves (Figs. 2 (a), 2 (b)). When several candidates 
are selected in the same event, only the candidate with 
the lowest \AE— A_E s hift| value is considered ("best can- 
didate"). From Monte Carlo studies, this algorithm is 
found to give the best reconstruction efficiency and the 
lowest cross-feed rate between the different D^D^K 
modes; it is found to introduce no bias on the signal 
extraction, since the latter is performed from the toes 
spectra only. However, in Fig. 2, to avoid the bias on AE 
inherent to this method, AE spectra are shown without 
applying this selection. 

IV. EVIDENCE FOR B -*■ D ( * } D ( * } K 

The toes and AE spectra of the selected events are 
shown in Fig. 2 for the sum of all the decay modes, sep- 
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FIG. 2: The AE and toes spectra (a,c,e) for the sum of all the B° —* D^*'D^*'K modes and (b,d,f) for the sum of all the 
B + _^ d^d^K modes. (a,b): AE for 5.27 < toes < 5.29 GeV/ c 2 . (c,d): toes for | AE - AE sM it\ < 2.5a A E- (e,f): toes for 
AE > 50MeV (background control region). The curves superimposed on the toes spectra correspond to the background fits 
described in the text and the shaded regions represent the background in the signal region 5.27 < toes < 5.29 GeV/c 2 . 



arately for B° and B + . The AE spectra are shown for 
events in the signal region defined by 5.27 < toes < 
5.29 GeV/c 2 . Signal events appear in the peak near MeV 
when reconstructed correctly, while the peak around 
-160 MeV is_due to D*DK and DD*K decays recon- 
structed as DDK and to D*D*K decays reconstructed 
as D*DK or DD*K. The toes spectra for the signal 
region arc shown for events with AE within ±2.5cta.e of 
the central AE value for the signal. The resolution ctab 
is determined from the data and is equal to 9.9 MeV for 
events involving no D*° and 11.3 MeV for events involv- 
ing one D*°. For events with two D*° candidates, the 
resolution is estimated from the Monte Carlo simulation 
to be 13.8 MeV. As explained above, only the candidate 



with the lowest | AE — A_B s hift| appears in the Toes spec- 
tra in case of multiple candidates. Both the toes spectra 
for the AE signal region and the AE spectra show clear 
evidence of a signal. On the contrary, the toes spectra 
for the background control region AE > 50 MeV do not 
contain any excess of events in the B signal region as 
expected. When fitting the Toes spectra, the combina- 
torial background component is empirically described by 
a threshold function [11] (henceforth referred to as the 
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ARGUS distribution), 



dN 
dm ES 



= /(m E s; AO = Am E s\h- - 



rn 



ES 



(4) 



x cxp 



-ci 



'ES 



where too represents the kinematic upper limit and is 
held fixed at the center-of-mass beam energy E^ eam = 
5.291 GeV, and A is a normalization factor. The func- 
tion depends on a free parameter £ that is determined 
from a fit to the toes spectrum of the background con- 
trol region. The number of combinatorial background 
events in the signal region is then estimated by nor- 
malizing the ARGUS distribution to the region 5.22 < 
Toes < 5.27 GeV/ c 2 in the AE slice containing the sig- 
nal (Figs. 2(c), 2(d)) and extrapolating it to the signal 
region 5.27 < to es < 5.29GeV/c 2 . The fitted ARGUS 
distributions are overlaid on the toes spectra of Fig. 2. 

The number of background events predicted in the sig- 
nal region by the fit is 1889±24 for neutral B mesons 
and 2512±27 for charged B mesons, while 2712 and 3482 
events are observed, giving an excess of 823±57 B° and 
970±65 B + events in the signal region. 



number of background events, /i° k s , in the signal region 
for this submode is calculated as 



. bk g 



5.29 



/ (x;A lk Xik)dx. 



(5) 



5.27 



If n k submodes arc used for a given mode, the branching 
fraction for that mode is then extracted by maximizing 
the following likelihood: 



n 

i=l 



N ik \ 



(6) 



where Nn~ and [ii k are the observed and predicted num- 
ber of events, respectively, in the signal region. [i ik is the 
sum of three contributions: 

• the predicted signal nf k , which is the product of the 
(unknown) branching fraction B k of decay mode k, 
the reconstruction efficiency e ikl the intermediate 
branching fractions Bf )D , and the number of BB 
events, N BB , assuming that the number of B°B° 
meson pairs produced at the T(4S) resonance is 
equal to the number of B + B~ pairs: 



fif k =B k xN BB xe ik x B\ 



DD. 



(7) 



V. DETERMINATION OF BRANCHING 
FRACTIONS 

In the following, the subscript k will be used to iden- 
tify the different B — > D^D^K decay modes (i.e., 
D°D°K+, D*-D°K+, ...). The subscript i will be used 
to identify the different decay submodes of the DD pair 
(i.e., i —Ktt x Ktt, Ktt x Ktttt , Ktt x K3n, ...). The 
subscript ik will therefore refer to B mode k decaying 
into DD submode i. 

The toes spectra obtained after a ±2.5<jae selection 
on (AE - AE shiit ) for all the different D^D^K modes 
are shown in Fig. 3 (B° decay modes) and Fig. 4 (B + de- 
cay modes). The corresponding event yields, computed 
as explained below, are given in Table I. In Figs. 3 and 
4 and in Table I, for a given B decay mode the sig- 
nals from the different DD decay submodes have been 
summed. However, to take advantage of the different 
signal-to-background ratios of the various submodes, the 
information from each submode is entered separately in a 
likelihood function used to calculate the B — > D^D^K 
branching fractions. As a first step, the ARGUS distri- 
bution shape parameter of each submode, Q k , is deter- 
mined from a maximum likelihood fit to the toes spectra 
of the background control region AE > 50MeV. An 
ARGUS distribution with the shape parameter £ fixed 
to this value Q k is then fitted to the toes distribution 
for the signal region \ AE — A£? s hift| < 2.5cta_e, excluding 
from the fit events with 5.27 < to es < 5.29GcV/c 2 . The 
factor A ik is calculated so that the function is normal- 
ized to the total number of background events and the 



• the number of combinatorial background events, 
Atj fe s , determined as described above (Eq. 5); 

• the peaking background /J?? 3, from other B — > 
£)(*)£)(*) x decay modes. The cross-feed between 
different DD decay submodes is found to be negli- 
gible and /^ ak is therefore calculated as 

/C* = x N BB x e'iil - ik) x B? D , (8) 

where e'(il — ► ik) is the cross- feed matrix ele- 
ment that represents the probability for B mode 
/ to be reconstructed as B mode k for DD decay 
submode i. The only significant cross- feed is ob- 
served between decay modes where a fake D*° re- 
places a real D* + or a real D°, for instance between 
D*-D°K+ and D*°D°K+, or between D*°D°K+ 
and D°D*°K+. 

The branching fractions B k for the sets of decay modes 
that have significant cross-feed are simultaneously fit, by 
maximizing the product FJ fe L k of the corresponding like- 
lihood functions. 

The D* and D branching fractions used in the branch- 
ing fraction calculation are summarized in Table II [12]. 
Branching fractions for decay modes reconstructed with 
a are calculated for neutral K mesons, including K° L . 
The selection efficiencies and the cross-feed matrices for 
each mode are obtained from a detailed Monte Carlo sim- 
ulation, in which the detector response is modeled with 
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FIG. 3: The tties spectra of the ten B° — > D^D^K modes. For each mode, all the D decay submodes used in the analysis 
have been summed, except for B modes for which the D x D decay mode is listed explicitly on the plot. The curves correspond 
to the background fits described in the text and the shaded regions represent the background in the signal region. Upper 
left: pure external W^-emission (spectator) decays B° — > D^*'° K + . Upper right: external+internal W-emission decays 

B° — > I)(*) - Z)(*)+ Kg. Lower left: pure internal VU-emission (color-suppressed) decays B° — > K s . 



the Geant4 program [13]. The simulated event samples 
of B — ► D^D^K decays used for the efficiency calcu- 
lation are generated according to a phase space model. 
For each decay submode, samples of about 15000 signal 



events have been produced. In addition, data are used 
whenever possible to determine detector performance: 
tracking efficiencies are determined by identifying tracks 
in the silicon vertex detector and measuring the frac- 
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FIG. 4: The 771es spectra of the twelve B + — > D^*' D*K modes. For each mode, all the D decay submodes used in the analysis 
have been summed, except for B modes for which the D x D decay mode is listed explicitly on the plot. The curves correspond 
to the background fits described in the text and the shaded regions represent the background in the signal region. Upper 
left: pure external W^-emission (spectator) decays B + — + K$ . Upper right: external+internal W-emission decays 

B+ -> d(*)°d ( -* )0 K+. Lower left: pure internal VK-emission (color-suppressed) decays B + — > D^*' D^*' + K + . 



tion that is well reconstructed in the drift chamber; the 
kaon identification efficiency is estimated from a sample 
of D*+ -► D°tt+, D° -> K-TT+ decays; the 7 and tt° ef- 
ficiencies are measured by comparing the ratio of events, 



N(t+ -► V T h+Tr°)/N(T+ -► I7 T /i+7r 7r ), to the published 
branching fractions [14]. Typical efficiencies range from 
20% for B+ -> D a D°K+ with both D° mesons decay- 
ing to K-TT+, to less than 1% for B+ ->■ D*~ D*+ K+ 
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TABLE I: Number of events and branching fractions for each mode. The excess (column 4) is the difference between the total 
yield observed in the signal region 5.27 < thes < 5.29 GeV/c 2 and the combinatorial background. It includes the contribution 
from the signal itself and from the cross-feed from the other D^D^K modes. The number of cross-feed events, computed 
from the cross-feed matrix and from the measured B — > D^D^K branching fractions, is given in column 5. When omitted, 
the predicted number of cross-feed events is smaller than 5 and has been neglected in the branching fraction calculations, 
as explained in the text. The first error on each branching fraction is the statistical uncertainty and the second one is the 
systematic uncertainty. For the decay modes with a significance S/VB smaller than 4, a 90% confidence level (C.L.) upper limit 
is also derived. Here, B is the sum of the combinatorial background and of cross-feed, while S = N — B, where N is the total 
yield in the signal region. The decay mode B + — > D*°D°K + has a large number of signal events but its significance is lower 
than 4 because of the large cross- feed from B + -» D°D*°K + and B° -» D*~ D° K + . The fractional statistical error on the 
branching fractions cannot be directly related to the fractional statistical error on the excess since the different decay submodes 
of the DD pair (not detailed in the table) enter with different statistical weights in the branching fraction calculation, while the 
yields given here are a raw sum over all the DD decay submodes. The statistical uncertainty on the background, dominated 
by the uncertainty on the ARGUS shape parameter is incorporated in the systematic error on the branching fractions. 



B decay 
mode 



Total yield 

N in the 
signal region 



Estimated 
combinatorial 
background 



Excess = 
Signal + 
cross-feed 



cross- 
feed 



Branching 
fraction 
(%) 



90% C.L. 

upper 
limit (%) 



B" -> D-D U K+ 
B" -> D~D*°K + 
B° -> D*~D°K+ 
B° -> D*-D*°K+ 



B decays through external W-emission amplitudes 



599 


479±12 


120±27 




0.17 ±0.03 ±0.03 


468 


337±10 


131±24 




0.46 ±0.07 ±0.07 


584 


399±11 


185±27 




0-3lloo3 ± 0.04 


289 


84±5 


205±18 




1.18 ±0.10 ±0.17 



B decays through external±internal W-emission amplitudes 





-> D~D+K U 


26 


19±2 


7±5 




0.08+005 ± °- 03 


0.17 


B° 


-> D*~D + K° ± D~D*+K° 


84 


34±3 


50±10 




0.65 ±0.12 ±0.10 




B° 


-> D*-D* + K° 


116 


48±4 


68±11 




0.88±g H ± 0.13 




B u decays through internal W-emission amplitudes 


B u 


-> D U D U K U 


175 


173±7 


2±15 




0.08 ±0.04 ±0.02 


0.14 


B° 


D°D*°K + D*°D°K 


248 


225±8 


23±18 




0.171 o^g ±0.07 


0.37 


B° 


-> D*°D*°K° 


123 


81±6 


42±13 


19.8 


0.33+^ ±0.14 


0.66 


B + decays through external W-emission amplitudes 


B+ 


D°D+K° 


367 


317±9 


50±21 




0.18 ±0.07 ±0.04 


0.28 


B+ 


- D*°D + K t] 


216 


175±7 


41±16 


9.6 


0.4ltg;^ ± 0.08 


0.61 


B+ 


D°D*+K° 


77 


31±3 


46±9 




0.52l° : ^ ± 0.07 




B+ 


-> D^'D'+K 


89 


43±4 


46±10 


9.0 


0.78^ ±0.14 





B + decays through external±internal W-emission amplitudes 



B+ 




D U D"K+ 


627 


469±11 


158±27 




0.19 ±0.03 ±0.03 




B+ 




D*°D°K+ 


552 


411±11 


141±26 


75.3 


0-18±8;8e ±0 - 04 


0.38 


B+ 




D°D*°K+ 


623 


402±11 


221±27 


37.1 


0.47 ±0.07 ±0.07 




B+ 




D*°D*°K + 


675 


468±15 


207±30 


66.6 


0.53l!5;iJ ±0.12 








B + decays throug 


i internal W-emission amplitudes 




B+ 




D~D + K + 


64 


65±4 


-1±9 




0.00 ±0.03 ±0.01 


0.04 


B+ 




D-D* + K+ 


45 


39±4 


6±8 




0.02 ±0.02 ±0.01 


0.07 


B+ 




D*-D+K+ 


64 


32±3 


32±9 




0.15 ±0.03 ±0.02 




B+ 




D *- D *+ K + 


83 


60±4 


23±10 




0.09 ±0.04 ±0.02 


0.18 



(D*+ -> D°ir+, D*- -» D°tt- ) with D° mesons decay- 
ing to K ~7r + 7T or i^~7r+7T~7r + . 



VI. SYSTEMATIC STUDIES 

Due to the large number of K ± mesons and to the 
large track multiplicities involved in the decays B — > 
£)(*)£)(*) iif, the dominant systematic uncertainties come 
from our level of understanding of the charged kaon iden- 
tification and of the charged-particle tracking efficien- 



cies. Both systematic uncertainties are estimated for 
each track and are given in Table III. Another impor- 
tant systematic is the uncertainty linked to the back- 
ground description. One of its components is from the 
uncertainty on the number of background events and is 
dominated by the uncertainty on the ARGUS shape pa- 
rameter £. The relative error on the branching fractions 
associated with this component varies from 5% up to 
20% depending on the mode and is uncorrelated from 
one mode to another. The other component is from the 
end point toq of the ARGUS distribution. Changing mo 
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TABLE II: Submode branching fractions used in the anal- 
ysis [12]. The errors on B(D° -> R-tt+tv ) and B(D° -> 
K~Tr + Tr~TT + ) correlated to the error on B(D° — * K~tt + ) are 
indicated separately with the subscript Kn. 



Mode 


B (%) 


£>" - K-n+ 

D° -> K-7V + 7V 


3.80 ±0.09 
13.10 ±0.84 ±0.31/^ 
7.46 ± 0.30 ±0.18k^ 


D+ -> K-n + iT + 


9.1 ±0.6 


D* + -> D+7T 


67.7 ±0.5 
30.7 ±0.5 


£*° - D° 7 


61.9 ±2.9 
38.1 ± 2.9 


-> n + n- 


68.60 ±0.27 



by ±1 MeV results in a ±1.4% change of the fitted com- 
binatorial background. The associated fractional error 
on each branching fractions is estimated to be ±3.5% in 
average and is correlated between all the modes. Other 
systematic uncertainties are due to uncertainties on the 
D and D* branching fractions, the 7r° reconstruction ef- 
ficiencies, the D vertex fit quality requirements, and the 
AE resolution used to define the signal box, as well as the 
statistical uncertainty on the efficiency due to the finite 
size of the Monte Carlo simulation samples and the un- 
certainty on the number of BB events in the data sample. 
The different contributions to the systematic uncertain- 
ties on the branching fractions are summarized in Table 
III. 

Possible decay model dependences of the efficien- 
cies were also studied by generating the decays B a —> 
D*~D+ and B° -» D*~D'+ (D+,D'+ -» D*°K+), 
where is the narrow (r = 1 MeV, m — 

2535.35 MeV/c 2 ) orbitally excited 1+ state of the D sJ sys- 
tem and D'+ is a wide (T = 250 MeV, m = 2560MeV/c 2 ) 
D s j resonance. The efficiency for reconstructing these 
modes was compared to the efficiency found for B Q — > 
D*~ D*°K + decays generated with a phase-space model. 
We found no statistically significant difference in efficien- 
cies; we assign a systematic uncertainty equal to the sta- 
tistical error of the ratio (5%). 



VII. SEARCH FOR RESONANT 
SUBSTRUCTURE 

B — > £)(*) £)(*)" K + decay modes are used to probe the 
possible presence of intermediate D s j resonances decay- 
ing into D^°K + , where D s j are P-wave excitations of 
the cs system. In the heavy-quark (charm) mass limit, 
the spin of the heavy quark decouples, and both the spin 
J of the meson and the total angular momentum (spin 
plus orbital) j q of the light quark become good quan- 
tum numbers [15, 16]. There are four P-wave states 
with the following spin-parity and light-quark angular 



momenta: 0+ (j q = 1/2), 1+ (j q = 1/2), 1+ (j q = 3/2), 
2 + (j q = 3/2). The two j q = 3/2 states can only un- 
dergo D-wave decay, and therefore have narrow widths. 
The remaining j q = 1/2 states decay via S- waves and are 
expected to be quite broad. Their masses are predicted 
to be ss 2.48GeV/c 2 (0+) and « 2.55GeV/c 2 (1+), while 
their widths are predicted to be a few hundred MeV [17]. 
However, the recent observation by the BABAR collabora- 
tion of a narrow state decaying to Dfir , with a mass of 
(2316.8 ± 0.4) MeV/c 2 (statistical error only) [19], would 
contradict these predictions and could indicate that the 
J p = + state has a mass lower than the D^K thresh- 
old; if this interpretation is confirmed, the + state would 
therefore not contribute to the B — > D^D^K final 
state. 

In the analysis described below, the two narrow reso- 
nances. D+ 1 (2536) and £>+ 7 (2573), are considered. The 
full Dalitz plot for the decay B° -> D*-D*°K+ is also 
examined. 



A. £>+(2536) 

(2536) is the most probable resonance to contribute 
to B — > D^D^K decays. It has already been ob- 
served and its measured parameters are m = (2535.35 ± 
0.60)MeV/c 2 , T < 2.3 MeV, J p = 1+ and j q = § [12]. 
Because of conservation of parity and angular momen- 
tum, only the decays £^(2536) — > D* K are allowed. 
In this analysis, a search is made for the 1)^(2536) in 
the final state D*°K + in the four decay modes B° — > 
D-D*°K+, B° -> D*-D*°K+, B+ -> 'd°D*°K+, and 
B + — ► D*°D*°K + . This resonance is not reconstructed 
in the D* + K® final state due to its lower reconstruction 
efficiency. 

Figure 5(a) shows the distribution of the variable 
Am = m (D*°K + ) - m (D*°) for the events recon- 
structed in the signal region (5.27 < rn.ES < 5.29 GeV/c 2 ) 
for these four decay modes. The distribution is fitted 
with a Gaussian function describing the signal. The com- 
binatorial background is represented by a threshold func- 
tion defined as 

g (Am) = N (Am - Am f e «( A ™" A ™o) . (9) 

The parameters of the Gaussian function (mean value 
Ami = 527.9 MeV/c 2 and standard deviation a Am = 
3.5MeV/c 2 ) are fixed to the values obtained from a fit 
to the same distribution resulting from the reconstruc- 
tion of inclusive .0^(2536) — > D*°K + decays in a large 
sample of events. This procedure yields an estimated sig- 
nal of 28t^ L»+ 1 (2536) -» D*°K+ events out of 764 ± 50 
B -» D^D*°K+ events. 

In order to extract upper limits on the contribution of 
£>+(2536) toB^ D^D^K decays, the same method 
is applied to the four individual decay modes, as shown 
in Fig. 6(a). The region 519 < Am < 537MeV/c 2 , il- 
lustrated by the dashed lines in Fig. 5(a), is defined as 
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TABLE III: Fractional systematic uncertainties on efficiencies and branching fractions. 



Item 


Fractional uncertainty on efficiency or branching fraction 


Charged track reconstruction 


0.8% per track for tracks with more that 12 hits required in the Drift Chamber 
1.2% per track for tracks without Drift Chamber requirement 


reconstruction 


2.5% per K", added in quadrature to the track reconstruction error 


7r^ reconstruction 


5 1% Der 7r ,J 


7 from ZT U -> D u 7 


5.1% per 7 (correlated with the 7r u systematic) 


1\ ILlUIlllllL-dLllJIl 


9 ^% npr K ± 


Vprt p"v 

reconstruction 


1 ^% npr *? trfipk" vprfp"v 

_1_ ■ O / U ul cIVj-LV V CI LCA 

3.1% per 3 track vertex 
5.7% per 4 track vertex 




2% for modes with zero or one D* u 
5% for modes with two D*°'s 


Background description 


5% to 20% (ARGUS shape parameter mode dependent) 
3.5% (end point mo) 


Monte Carlo statistics 


2% to 10% per DD submode (mode and submode dependent) 


Intermediate br. fraction 


see Table II 


Number of BB 


1.1% 


Decay model 


5% 



> 

s 




0.5 0.55 0.6 0.65 0.7 0.75 



m(D , * ) V)-m(D , * )0 ) (GeVAf) 



FIG. 5: (a) Am = m (D*°K+) - m (D*°) distribution for 
events reconstructed in the B — > D 1 -*' D*° K + signal regions, 
(b) Am = m (D°K + ) —m (-D ) distribution for events recon- 
structed in the B —* D^*' D° K + signal regions. 



suits obtained and gives a 90 % confidence level (C.L.) 
upper limit on the quantity i?2536: 



-R2536 



n{B -> D«£>+(2536)) x B (£>+ (2536) -> D*°K+) 

n (B -> D(*W*°K+) ' 

(10) 

where n(B -» D W .D+ (2536)) x B (£)+ (2536) -> D*°K+) 
and n(B — » D^D*°K + ) are the observed num- 
ber of signal events. Using the B — > _D(*)Z?* if + 
branching fraction measurements of Table I, these re- 
sults can be compared to the only existing measure- 
ment of inclusive 1)^(2536) production in B decays, 
B (B -> £>+(2536)X) < 0.95% at 90% C.L. [18]. 



TABLE IV: ££(2536) 
D M D*°K+ decays. 



D*°K + contributions to B 





Total yield 




£^2536 


B decay mode 


in £>+(2536) 


Estimated 


90% C.L. 




signal region 


background 


upper limit 


D-D* U K+ 


16 


7.8 ±0.6 


11.6% 


D*~D*°K + 


13 


7.3 ±0.6 


5.8% 


D°D*°K+ 


12 


11.1 ±0.8 


3.1% 


D*°D*°K+ 


20 


8.7 ±0.5 


9.1% 



the signal region and the number of combinatorial back- 
ground events in this region is estimated from the fit by 
integrating the background function g defined in Eq. 9. 
The total number of events observed in the signal box 
is compared to the expected combinatorial background 
when extracting the limits. Table IV summarizes the re- 



B. 



D+ (2573) 



The contribution of the D+j(2573) resonance to B — ► 
D(*)D { *~>K decays is also studied. This resonance is 
thought to be the other narrow state in the j q = 3/2 or- 
bitally excited D s j doublet, together with the ££(2536). 
The world average values of its mass and width are 
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FIG. 6: (a) Am = m (D*°K + ) — m (D*°) distributions of the events reconstructed in the four B — > D^*' D*° K + decay 
modes, with 5.27 < ities < 5.29GeV/c 2 . (b) Am = m (D 7f + ) — m (-D ) distributions of the events reconstructed in the four 
B — > D^*> D°K + decay modes, with 5.27 < wes < 5.29 GeV/c 2 . These distributions are fit with the sum of a threshold function 
/ (Eq. 9) for the background components and a Gaussian function for the D s i(2536) and D s j(2573) components. The mean 
value and standard deviation of the Gaussian distributions have been fixed to the values obtained from a fit to the inclusive 
D s i(2536) and D s j(2573) samples, as described in the text. 



m = (2573.5 ± 1.7) MeV/c 2 and T = (15±|) MeV [12]. Its 
spin-parity has not been measured but its natural width 
and decay properties are consistent with a J p — 2 + state 
[12]. If it is indeed a spin-2 resonance, it cannot be ob- 
tained with a VF-mediated tree diagram but might still 
be reached through final state interactions. 

The allowed decay modes of the Dfj (2573) are DK 
and D*K, both proceeding through a D-wave. Because 
of the limited phase space, the latter is highly suppressed 
[17]. In this analysis, a search is made for the £)j"j(2573) 
in the decay mode D°K + , in the four channels B° — > 
D-D°K+, D*~D°K + , D°D°K + , and D*°D°K+. 

The method developed for the D^(2536) study is ap- 
plied. Figure 5 (b) shows the Am = m (D°K + )—m (D°) 
distribution for the events reconstructed in all four B — > 
D(*)D°K+ decay modes. The mean value and the stan- 
dard deviation of the Gaussian component of the fit func- 
tion are fixed, respectively, to Am = 708 MeV/ c 2 and 
CAm = 6MeV/c 2 , which are the values derived from a 
large inclusive D^j (2573) — > D Q K + data sample. The 
fitted yield of £>+ 1 (2573) -> D°K+ decays is 13±9 events 
out of 604 ± 54 B -> D^D°K + events. 

Defining the signal region 687 < Am < 729MeV/c 2 , 



90% C.L. upper limits on the contribution of D^j (2573) 

to B -> DWDWjr decays are set for each of the four in- 
dividual decay modes (Fig. 6(b)). The number of events 
observed in the signal box, the number of background 
events expected from the fits, and the resulting limits 
are given in Table V. -R2573 is defined here as 

n(B -> (2573)) x B (D+ (2573) -» D°K+) 

JT 12573 = 



B -> D(*W°K+) 



(11) 



TABLE V: _D+ 7 (2573) 
Di*) D o K + decays 



D°K + contributions to B 





Total yield 




^2573 


B decay mode 


in D+ (2573) 


Estimated 


90% C.L. 




signal region 


background 


upper limit 


D~D U K + 


25 


26 ±3 


6.3% 


D*~D°K + 


41 


42 ±3 


5.0% 


D°D°K+ 


38 


29 ±3 


12.2% 


D*°D°K + 


37 


30 ±3 


12.3% 



16 



C. Dalitz-plot analysis of the decay 
B° -». D*~D*°K + 

As suggested in Rcf. [20], the study of decays B — > 
£>(*)£)(*)K could be used to search for evidence of the 
yet undiscovered broad j q = 1/2 Z? s j states, if the de- 
cays D s j — > D^K are allowed by the available phase 
space. The decay mode £?° — > D*~ D*°K + , which has 
the largest number of reconstructed events and also has 
the largest purity, is used for this search. The results 
are shown in Fig. 7. The upper right plot (Fig. 7 (b)) is 
the Dalitz plot m 2 (D* Q K+) vs. m 2 (D*-D* a ) expected 
for three-body B° — > D*~D*°K + decays generated with 
a phase space model. The Dalitz plot m 2 (D* a K + ) 
vs. m 2 (D*~ D*°) for data events in the signal region 
5.27 < m ES (D*-D*°K + ) < 5.29GeV/c 2 (Fig. 7 (a)) is 
shown in Fig. 7 (c). The next three plots (Figs. 7 (d), 7 
(c), 7 (f)) show the projections m(D*°K+), m(D*-K + ), 
and m(D*~ D*°) for the same events. The hatched his- 
tograms show the contribution expected from the com- 
binatorial background; their shapes are derived from the 
events with m ES (D*- D*°K+) < 5.26GcV/c 2 . The open 
histograms show the contribution expected for three- 
body B° — > D*~D*°K + decays generated with a phase- 
space model. The density of events in the lower region 
of the Dalitz plot (i.e., for small values of m(D*°K + ) 
and large values of m(D*~ D* )) is significantly larger 
in the data ( Fig. 7 (c)) than in the simulation with 
no resonance ( Fig. 7 (b)). It could be interpreted as 
the presence of a broad resonance decaying to D*°K + , 
like the J p = 1 + , j q = 1/2 state predicted by Heavy 
Quark Symmetry models [15-17]. However, more events 
are necessary to confirm this hypothesis and to estimate 
the resonance properties such as mass and width. 

As previously discussed (Sec. VI), the hypothetical 
presence of broad resonances in the decay chain is ac- 
counted for by a 5% relative systematic error on all the 
B — > D^D^K branching fraction measurements de- 
scribed in this paper. 

VIII. CONCLUSIONS 

A measurement of the branching fractions for the 22 
B -> D^D^K modes is given in Table I. For the 
decay modes for which S/VB is smaller than 4, a 90% 
C.L. upper limit is also derived (here, B is the sum of the 
combinatorial background and the cross-feed background 
from other D^D^K modes and S — N - B, where N 
is the total yield in the signal region). This is the first 
complete measurement of all possible B — > D^D^K 
channels. The measured branching fractions are in good 
agreement with earlier measurements made with smaller 
data sets for some of these modes [5-8] . 

The existence of the decays B a — > D*~D + K° S and 
B a — > D*~ D* + K°, which are an admixture of CP even 
and CP odd eigenstates, has been demonstrated. These 



decay modes could be used in the future, with larger 
event samples, to also determine sin2/3 and cos2/3 [21- 
23]. A significant signal for the color suppressed decay 
mode B + — > D*~D + K + has also been observed. 

One of the motivations of this analysis is to understand 
whether decays B -> D^D^K can explain the wrong- 
sign D-meson rates in B decays and reconcile the total 
b — > ccs rate with the predictions of Ref . [3] . After sum- 
ming over all submodes, the branching fractions of the 
B° and of the B+ to £)(*)£>(*) A" arc found to be 

B(B° -> D M D M K) = (4.3 ± 0.3(stat) ± 0.6(syst)) %, (12) 

B(B+ D^D^K) = (3.5 ± 0.3(stat) ± 0.5(syst)) %. 

(13) 

This study shows that a significant fraction of the 
transitions b — > ccs proceed through the decays B — > 
£)(*)£)(*) A. These decay modes account for about one 
half of the wrong-sign D production rate in B decays, 
B(B -» DX) = (7.9 ±2.2)% [4]; however, because of the 
large statistical error on the latter measurement, it is not 
yet clear whether they saturate it. 

A search for resonant substructures shows that the 
££(2536) contribution to B -» D^D*°K + decays is 
small. No evidence for a ££(2573) contribution to 
B — ► D^D°K + decays is found. Finally, a simple 
Dalitz-plot analysis of the decays B° -> D*-D*°K+ 
shows that the three-body phase-space decay model does 
not give a satisfactory description of these decays. 
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FIG. 7: Dalitz plots and projections for the decay B° — ► D* D*°K + . The content of the different plots is discussed in the 
text. 
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